Nonlinear dissipation can combat linear loss 
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We demonstrate that it is possible to compensate for effects of strong linear loss when generating 
non-classical states by engineered nonlinear dissipation. We show that it is always possible to 
construct such a loss-resistant dissipative gadget in which, for a certain class of initial states, the 
desired non-classical pure state can be attained within a particular time interval with an arbitrary 
precision. Further we demonstrate that an arbitrarily large linear loss can still be compensated by 
a sufficiently strong coherent or even thermal driving, thus attaining a strongly non-classical (in 
particular, sub-Poissonian) stationary mixed states. 



PACS numbers: 03.65.Yz, 42.50.Dv 



Nowadays, an engineered dissipation for quantum state 
manipulation is an intensely developing field. More that 
a decade ago, it has been shown that in the systems such 
as ions in magnetic traps it is possible to tailor non- 
linear dissipation in a rather wide range [1], and later 
the concept of the "quantum state protection" was born 
0. Different kinds of a nonlinear dissipative apparatas 
(aptly nicknamed "dissipative gadgets" @, 3) have been 
shown to be useful for many important tasks, for ex- 
ample, for generating non-classical and entangled states 
of few-body and many-body systems (5T-fr2l]. performing 
universal quantum computation 0, [HJ , constructing dis- 
sipatively protected quantum memory performing 
precisely timed sequential operations, conditional mea- 
surements or error correction jl 31 ] - The central idea of 
all the dissipative gadgets is to make dissipation to drive 
the system towards a desired steady state (which is prac- 
tically independent of the initial state). Remarkably, in 
this context, dissipation serves as a helpful quantum re- 
source rather than being a hindrance. 

However, in all these currently known schemes, the 
engineered dissipation is far from being an universal effi- 
cient tool for combating the usual linear loss, inevitably 
present in any realistic dissipative gadget. For example, 
the conventional single-photon loss makes the generation 
of a non-classical pure stationary state impossible for just 
any kind of nonlinear dissipation [3]. Just as it is for 
the coherent control, one is generally obliged to mini- 
mize linear losses using some extra effort while arranging 
for the nonlinear terms to produce a desired effect (see, 
for example, Ref. || In particular, it is rather 

hard to produce a desired state of an electromagnetic 
field for schemes relying on optical nonlinearities 043- 
However, surprisingly enough, the nonlinear dissipation, 
when properly designed, do can combat the effects of an 
arbitrarily strong linear loss, for finite time intervals and 
for stationary states. That is the main message of our 
contribution. 

We start with demonstrating that it is always possi- 
ble to generate a state approximating the desired pure 
non-classical state with any given precision for an arbi- 
trary ratio of linear and nonlinear loss rates. To illustrate 



our argument, let us consider a simple Lindblad master 
equation for the single mode of the electromagnetic field 

4/?= -i[H,p]+T(n+l)C{a)p + TnC{a f )+-/C{A)p (1) 
at 

where the nonlinear dissipation is described by the Lind- 
blad operator A, the operators a', a are the usual bosonic 
creation and annihilation operators, and the operator H 
represents the system Hamiltonian. T and 7 > are rates 
of linear and nonlinear dissipation, and the parameter n 
represents the average number of photons in the thermal 
pump. The Liouvillian C(x) acts on the density matrix 
as C(x)p = 2xpx> — x*xp — px^x. 

It has been proved that the coherent state is the only 
possible pure stationary state of Eq.JT]) |3| for T > 
(and it will be the vacuum state if no coherent driving 
is present). Still there is a possibility to generate any 
desired state during the time interval t <C To ilus- 
trate the principle, consider A = \(f>) (y\a k , k > 1, as the 
Lindblad operator for the engineered dissipation, where 
the vector \y) describes the initial state and \4>) is the 
target pure state. The density matrix satisfying Eq.(p} 
with the above nonlinear A-term has the property: 



Pp(t)P = Tr{Fp(0)}|0)(0| + 0(e 



(2) 



Here 7 off = min{l,2(l - \(^)\ 2 )}Afj with *) = 
(N)~ 1 / 2 (a)) k \y). The norm-factor is given by Af = 
{y\a k {a^) k \y). The state |0) = (1 - K0|*}| 2 ) -1/2 (|$ - 
(\I r |^>)|\I r )) is the orthogonal complement of \*f>) in the sub- 
space spanned by \<p) and \^S), with P being the projector 
on this subspace. In particular, Eq.@ gives 



p^M = (l-|(0|vI/)| 2 )Tr{P jO (O)}, 



(3) 



that is, for Pp(0) = p(0) the large Af in Eq.© im- 
plies that the system will rapidly evolve to the state \<j>) 
with the fidelity 1 — K^l^)! 2 . Furthermore, note that if 
((fil^) = 0, one can drive the system from an initial state 
to the state \<j>) orthogonal to it. 
For a linear loss present, a simple physical picture be- 
hind the loss-suppression mechanism is provided by sys- 
tem jumps to the lower energy levels. Notably, the rate 
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FIG. 1: (Color online) Generation of nonclassical states from 
a coherent state \a) in the presence of linear loss using de- 
signed nonlinear dissipation described by the Lindblad op- 
erator A. Thermal driving is absent, n — 0, and the loss 
ratio is given by 7 = T/5. I. Generation of the Fock state 
1 2) using nonlinear dissipation given by A — 2){a|a 2 . (a) 
Fidelity of the target two-photon state as a function of non- 
linear loss Ft. Solid, dotted, dashed and dash-dotted lines 
correspond to a — 2, 3, 4, 5. The panel (b) shows the photon 
number distribution of the generated state when the fidelity 
in (a) is maximal; black, dark grey, grey and light bars corre- 
spond to a = 2, 3, 4, 5. II. Generation of sub-Poissonian light 
and photon-number squeezing using tailored dissipation with 
A = a(a + a — 1). (c) Mandel Q-parameter vs nonlinear loss 
Ft; solid, dotted, dashed and dash-dotted lines correspond to 
a = 2,4, 6, 8. (d) Photon number distribution of the gen- 
erated state for the minimal Q and a = 8 (grey bars) in 
comparison with the Poissonian distribution with the same 
average number of photons (solid line). 

of transition to the lower levels due to nonlinear dissi- 
pation can be much higher than the transition rate due 
to linear loss. Now, the stationary state of the evolution 
induced by nonlinear dissipation corresponds to a non- 
classical state. Hence when the dynamics due to non- 
linear dissipation prevails over that due to linear dissi- 
pation, the system is driven into a desired nonclassical 
state with influence of linear loss being negligible. De- 
signing dissipation (i.e., the Lindblad operator) we deter- 
mine the form of the nonclassical target state. If we apply 
this quantum jump formalism to the system desribed by 
Eqs. (HJ)-©, we can deduce that for 

T(n + l)(*|a t a|*) < 7 (*|A t A|*) (4) 

the influence of the linear loss (and of incoherent driv- 
ing) during the time of the target state generation is 
negligible. For a coherent initial state \y) = \a) with 
sufficiently large amplitude, a 3> 1, condition Q gives 
r(n + 1) <C 7A/*. In practical terms, for the target state 
lying within the subspace of Fock vectors with up to k 
photons, the generation can be nearly perfect and im- 
proved further just by increasing the amplitude a. 

Fig. [T] (a,b) illustrates the solution of Eq.([T]) for the 
generation of the two-photon Fock state from an ini- 
tial coherent state |ev) with the nonlinear dissipation de- 
scribed by the Lindblad operator A = |2)(a|a 2 . The 



rate of nonlinear loss, V, is 5 times less than the rate 
of the linear loss, 7. The incoherent driving is absent. 
By increasing the amplitude a the target state \2) is ap- 
proximated with a larger fidelity and over a shorter time 
period, despite the presence of rather strong linear loss 
(FiglHa.b)). 

As we have seen, if for a class of states the jump rate 
of the nonlinear dissipation far exceeds the jump rate for 
the linear loss, then an influence of the linear loss is neg- 
ligible. Moreover, the dynamics due to the engineered 
nonlinear dissipation can also be non-exponentially fast. 
For example, for the two-photon dissipation, quite dif- 
ferent initial states can decay into the stationary state 
practically over the same period of time (16J. Thus, for 
an appropriately constructed dissipative gadget it is al- 
ways possible to combat the linear loss with any pre- 
defined rate just by choosing the appropriate initial state, 
in particular, by increasing an average number of pho- 
tons in the initial state. Of course, there arc diffircnt 
practical limitations in construction of dissipative gad- 
gets. However, even for available types of nonlinear dis- 
sipative processes it is possible to have significant non- 
classicality (in particular, large photon-number squeez- 
ing) on the time scales when influence of the linear losses 
is negligible (see, for example, Refs. d,[l7j])- ^his * s mus ~ 
trated in FigQJc, d), where the exact solution of Eq.(p} 
is shown for the nonlinear dissipation described by the 
Lindblad operator A = a(a + a — 1) in the presence of a 
linear loss and for an initial coherent state \a). Obvi- 
ously, this Lindblad operator satisfies condition (j4| for a 
sufficiently large a. In the initial stages of evolution, non- 
classical states can be generated and the non-classicality 
(photon-number squeezing in this case) is increased with 
increasing of the amplitude a in the initial state. Panel 
(c) shows the Mandel Q parameter confirming that the 
generated light exhibit the sub-Poissonian statistics (for 
Q < 0). Note, that the maximum of non-classicality is 
reached rather quickly. The generated state then further 
evolves towards the classical regime under the influence 
of the linear loss but it happens quite slowly compared 
to the time required for the nonclassical state generation. 
This is a rather general feature of the dynamics for any 
dissipative gadget satisfying condition §t§ fl8j . 

It might seem surprising, but analogous approach can 
allow for combating linear losses also in the long-time 
limit. More precisely, there is a class of nonlinear losses 
such that for an arbitrarily large (but finite) ratio e — 
r/7 an influence of losses on the stationary state can 
be completely eliminated by the sufficiently strong co- 
herent driving, and, moreover, even by thermal pump- 
ing. Note, that the steady state, while being mixed, 
can still be strongly non-classical. Consider, for the ex- 
ample, nonlinear coherent loss (NCL) described by the 
Lindblad operators A = af(a 1f a), f(x) being a smooth 
function. This operator is the annihilation operator for 
the so called "f-deformed" quantum harmonic oscillator; 
eigenstates of thi s op erator are reffered to as "nonlinear 
coherent states" Any pure state non-orthogonal to 
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FIG. 2: (Color online) Generation of nonclassical states in the 
long-time limit using nonlinear coherent loss (NCL) to coun- 
teract linear loss, (a) Mandel's Q-parameter vs amplitude of 
the coherent driving ao = f2/7 for the NCL with A — aa^a—l 
provided by the exact solution of Eq.fTJ; solid, dash-dotted 
and dashed curves correspond to the relative loss rate T = 7, 
57, IO7. Inset shows photon-number distribution of the sta- 
tionary state for ao — 150; light-grey, grey and black bars 
correspond to T = 7, 57, IO7. (b) Same for the NCL with 
A = a(a'a — 1) for the exact solution of Eq.fT]) (dashed line) 
and for the solution of the approximate equation ([6]) (solid 
line). Inset shows photon- number distributions for the solu- 
tion of the exact (grey bars) and of the approximate (black 
bars) equations, (c) Exact (dashed line) and approximate 
(solid line) solutions for Q-parameter for A = a(a^a — l) 2 ; the 
inset shows long-time behavior of Q-parameter given by the 
exact solution, (d) Exact solution for Q-factor in dependence 
of the average number of thermal pumping photons, «t; the 
coherent driving is absent; A = a(a^a — l) 3 . 



an arbitrary Fock states can be exactly represented as a 
nonlinear coherent state. If it is orthogonal to some Fock 
states, then one can still devise a nonlinear coherent state 
closely approximating the state in question [20| . NCL 
can be realized in practice in a number of schemes (for 
example, with ions or atoms in traps Bose-Einstein 
condensates 2l| or even in multicore nonlinear optical fi- 
bres @. With function f(x) having a finite or countable 
number of zeros, one can generate Fock states or even 
" comb" the initial state by filtering out some pre-defined 
set of components in the given basis [22j . 

Let us show now that even a weak NCL can be pro- 
tected from large linear loss and used for generation 
of non-classical stationary states. We use the classi- 
cal driving in the form H = ifl(a — cv) where the pa- 
rameter fi represents the strength of driving and, for 
simplicity, is taken to be real. For the moment be- 
ing, the thermal pumping is assumed to be absent, 
n = 0. Fig. [5] depicts such a non-classicality rescue 
procedure when generating sub-Poissionian light. No- 



tably, for A = a(a' i a — 1), the Mandel's Q parameter 
Q = (((a f a)) 2 ) - ((at a ))) 2 )/((at a ))) - 1 eventually tends 
to the same limiting value for quite different linear loss 
rates, r = 7, 57 and IO7 (Fig. [Ha)). It is remarkable 
that the generated states are also quite similar (see inset 
in Figl^a) for the density matrix). 

The nature of this phenomenon can be well illustrated 
and clarified with the help of the following simple ap- 
proximation. Eq.([T]) can be written as 



dt' 



jap(B^ - a ) + 7(.B - a )pa> 

ja?(B - a )p - 7/°(- Bt - a o) a ~ 
7a[[p,/(a t a)],/(a t a)]a t 



(5) 



where B = a([f(a^a)] 2 + e). The last term is small in 
comparison with the others when the inequality {f{n) — 
f(m)} 2 -C 2f(n)f(m) holds for all n and m, correspond- 
ing to essentially non-zero density matrix elements p nm - 
For example, for any power-law nonlinearity f(x) ~ x k , 
this condition is satisfied when the driving is sufficiently 
strong. Neglecting the last term in Eq.©, the station- 
ary state, p s , for Eq.© is the eigenstate of the opera- 
tor B satisfying Bp s = a n p s and p s B^ = a n p s ; where 
cio = This leads to the following recurrence rela- 

tion for the diagonal (in the Fock-state basis) elements 
of the steady state: 



Pn 



(6) 



Assuming that the steady-state photon number distribu- 
tion is peaked at no, from Eq.© we get a simple condi- 
tion for determining hq: 



Pn,n ~ Pn— l,n- 



n ((f(n ) 2 +e) 2 ) 



This condition indicates that for sufficiently strong 
classical driving and for function f(n) increasing mono- 
tonically for sufficiently large n > 0, one will always have 
/(no) 2 3> e, and the influence of linear losses on the 
generated steady state will be negligible. The relation 
© also allows for estimating the width of the photon- 
number distribution of the stationary state. Assuming 
that the value of the diagonal element p s n n changes only 
weakly for small changes Sn of the photon number n near 
no, one obtains from Eq.® the following expression: 



PriQ-\-5n,no+8n ^ Pr, 



exp 



+ ( , W(no)/(n ) \ 
2n I /(n ) 2 +£ j 



(7) 



From Eq.([7]) we can estimate the variance, A 2 n, of the 
steady-state photon number distribution: 



A 2 n « n 1 + 



4n /(no)/(n ) 
f(n ) 2 +e 



(8) 
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Thus, our approximation shows that in the limit of 
strong driving the steady state produced by the NCL 
will always be photon-number squeezed state, provided 
that /(no) > 0. The squeezing can be quite high. For 
example, let us use f(n) = n — 1, which has been proved 
feasible u sing three-well potential in Bose-Einstein con- 
densates (21| or three-core nonlinear fibers @. For the 
simple NCL with f(n) = n — 1 , the value A 2 n — > uq/5 
is asymptotically reached for n ~ ("o) 2 ^ 5 — > +oo. 
This value corresponds to the Mandel parameter equal 
to —0.8. Fig. [21b) shows that the approximation to 
Eq. (J5]) indeed gives rather good estimate for both the 
Mandel parameter and the generated state. It should 
be noted though, that for rapidly increasing /(n) the 
approximation works somewhat worse. It can be easily 
seen from Eq.©, since it predicts Q-parameter close to 
— 1 for rapidly increasing f(n), e.g., Q ~ — 4k/ (Ak + 1) 
for f(n) — n k . Q = — 1 corresponds to the Fock state. 
But when the photon number distribution becomes very 
narrow, the assumption of slowly changing p nn near max- 
imum of the photon number distribution can be hardly 
applied. Nevertheless, the approximation still provides a 
qualitatively correct description as illustrated in Fig.^c) 
for / = (x — l) 2 , where exact and approximate solutions 
for Q-parameter are compared. Hence even for f{x) in- 
creasing rather rapidly one can indeed have highly pro- 
nounced non-classicality unaffected by any linear loss. 
However, it happens at the expense of quite strong co- 
herent driving required (see Fig. [21c)). 

The essence of our method to preserve non-classicality 
lies in engineering the evolution to yield a stationary state 
obeying the following requirement. A stationary state of 
the master equation for a given system should be such, 
that the transition rates corresponding to nonlinear loss 
far exceed transition terms corresponding to linear loss, 
that is, the condition ((4]) holds. This simple fact points to 
the possibility of exploiting not only the coherent driving, 
but other kinds of driving, too. In particular, the conven- 
tional linear thermal driving can protect non-classicality 
as well, although this seems quite paradoxical. Indeed, in 
the absence of the coherent driving (f2 = 0), from exact 
Eq.© the following recurrence relation is obtained: 

n 

Pn = Pn - 1 (n + l) + ( 7 /T)P(n) (9) 

For small values of {"f/T)f 2 {n) the state © is very close 
to a thermal state. For (7 /T) f 2 (n) 3> (n + 1) the photon 



number distribution is effectively truncated. The trun- 
cation number, n t can be estimated from the following 
equation: n + 1 ~ (^f/T)f(n t ). If the density matrix 
elements p n decrease with increasing photon number n 
faster than for any coherent state: p n / Pn-i — o(l/n), 
the considered state is nonclassical. Eq.(j9]) implies that 
for f(x) growing faster than x this condition is satis- 
fied and the stationary state is nonclassical for any val- 
ues of system parameters. Fig. [2jd) shows an example 
of such a "thermal rescue" for the nonlinear loss with 
f(x) = (x — l) 3 . Thus, remarkably, the thermal excita- 
tion is able to produce photon anti-bunching. However, 
it is easy to see that the minimal Q-value always corre- 
sponds to some finite n. 

To conclude, we have shown that dissipative gadgets 
can be made extremely robust with respect to addi- 
tional linear losses unavoidably affecting any realistic en- 
gineered dissipation scheme. First we have shown that it 
is always possible to devise a dissipative gadget for the 
generation of any desired target state starting from the 
wide class of initial states during the time interval when 
influence of linear loss is negligible. We provided several 
examples illustrating such scenario using coherent input 
states. The larger the difference in average number of 
photons between the initial state and the target state, 
the more robust the scheme can be made. Then, by ap- 
plying conventional linear quasiclassical driving, the non- 
classicality of the stationary state generated by the dis- 
sipative gadget can be preserved. It is remarkable that 
there exists a class of dissipative gadgets based on the 
nonlinear coherent loss (NCL), for which the influence 
of linear loss can be completely annihilated by a suf- 
ficiently strong coherent driving. The resulted state is 
mixed. However, this stationary mixed state can be very 
close to the Fock state exhibiting strong photon-number 
squeezing. Furthermore, there exists a class of dissipative 
gadgets, for which non-classicality of the stationary state 
can be rescued even by using usual incoherent thermal 
driving. 
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